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Abstract 

The replicated field theory of the random field Ising model involves the 
couplings of replicas of different indices. The resulting correlation func- 
tions involve a superposition of different types of long distance behaviours. 
However the n = limit allows one to discuss the renormalization group 
properties in spite of this phenomenon. The attraction of pairs of replicas 
is enhanced under renormalization flow and no stable fixed point is found. 
Consequently an instability occurs in the paramagnetic region, before one 
reaches the Curie line, signalling the onset of replica symmetry breaking. 
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1 Introduction 



In a recent article we have briefly examined the field theory associated with the 
random field Ising model|l|, within the replica approach. It is a 4 field theory, 
but with the following modifications : 

(i) there are n fields <f> a , a — 1, • • • , n and the only symmetry is the permutation of 
the replicas. This allows for five couplings namely U\ V] </r, u 2 (fia^b, % 0^6; 

a* ab ab 

n 

UiJ2 ( t )2 aMc,U 5 (J2<Pa)' i - 
abc 1 

(ii) in the n = limit the (bare) propagator at T c is given by 

G ab { P ) = + (1) 

The second term results from the averaging over the random field h(x) : 

< h{x)h(y) >= A5{x - y) (2) 

At first sight, the l/p A singularity of the propagator could imply that the upper 
critical dimension is eight, instead of six, but it will be argued that it is indeed six 
because of the n = limit. Then the renormalization group studies which were 
conducted long ago || f|, [5[ dealt only with the single coupling constant ui, and one 
A/p 4 propagator per loop [contributions with less than one A/p 4 per loop being in 
any case infra-red subdominant]. This limitation was the result of considering the 
averaging over the random field of connected correlation functions. However once 
one introduces the coupling constants «2, • • • , «5, which involve several replicas, one 
must also consider contributions with moore tan one A/p 4 per loop. This modifies 
notably the conventional renormalization programme, and we shall carry it in some 
detail here in the zero replica limit. At the end we recover the n = limit of the beta- 
functions found in [|IJ and thus confirm its main conclusion, namely the instability of 
the dimensional reduction fixed point. 

Furthermore we are left with infra-red singular and attractive contributions to 
the 4-point function of type 113, i.e. corresponding to two distinct replicas. They can 
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give rise to negative eigenvalues in the iteration of the related Bethe-Salpeter kernel, 
implying the occurence of a glassy phase in the paramagnetic domain, before one 
reaches the Curie line ||. We shall briefly return to this point at the end. 

2 Upper critical dimension 

Let us first write the Boltzmann weight for this replicated field theory : 

/3H = J d d x{\ £ [(V0 a ) 2 + vMaf] ~ J E tafa 
a ab 

+^E^ + |E €h + 1 E m + x E + |(E **)*}. (3) 

i ab aft abc ' a, 

Indeed in [|IJ it was argued that the standard derivation of a field theory from the spin 
model in the presence of a random field, namely the consideration of the fluctuations 
around mean field theory, did yield those five coupling constants. The corresponding 
propagator at the critical temperature is thus simply given by 

Gab(p) = + ~jr^ tt)- (4) 

n-s-0 p z p z [p z — 71 A) 

In view of this A/p 4 one does get infrared singularities in dimensions lower than 
eight. For instance there is an obvious one- loop contribution to the renormalization 
of M3 proportional to (Aui) 2 which is singular at low external momenta as l/p 8 ~ d . 
However if one considers higher loops they are either less singular or they vanish 
with n. Indeed the diagrams which would be maximally singular in eight dimensions 
contain a A/p A on every internal line. Therefore as soon as the diagram contains any 
vertex which is not connected to one or two of the four external lines, which is bound 
to happen at most at order five in perturbation theory, it contains at least one free 
sum over replica indices, and thus it vanishes with n. 

Therefore the singularities that one encounters in these functions between eight 
and six dimensions are due to a finite number of graphs and thus there is no pos- 
sibility of anomalous dimensions, which could modify the behavior given by simple 
dimensional analysis. This is very much like, say a six-point function at criticality, in 

2 



a 4 theory : it is singular at low momentum for d < 6 but this singularity remains 
canonical down to four dimensions. Here similarly these singularities remain given 
by these few graphs down to d = 6, at which a full renormalization analysis becomes 
necessary : the upper critical dimension is six, not eight, because of the n = limit. 

3 Renormalization of the coupling constants 

Given the effective Hamiltonian (|3|), we wish now to renormalize the four-point 
functions, corresponding to the five coupling constants Uj(j = l,---,5), retaining 
terms involving one A-propagator per loop or more. We consider the one-particle- 
irreducible four point function F^} cd (pi, • • • ,Pa) and in order to minimize the number 
of momenta involved, we choose for simplicity the symmetric point 

Pi -Pj = -j (5<% - 2) (5) 

4 

compatible with momentum conservation = 0. We work with a dimension- 

i 

ally regularized theory in dimension d — 6 — e , and will renormalize by minimal 
subtraction. 

For simplicity let us first keep only u\ and alone and examine what is happen- 
ing. The four-point function of type 1 , namely the one which involves a product of 
three Kronecker deltas S a bS ac 5 a d, is given , at one-loop order, by 

r\ J (p) =m -3«i H . (6) 

e p e 

If we mutiply it by A, one sees that the real coupling constants which enter into all 
the diagrams which involve exactly one A per loop are 

9i = Auj. (7) 

Their (inverse length) dimension is e = 6 — d. Therefore in a minimal scheme, the 
relation between the bare gi and the dimensionless renormalized gf is, at this order 

5l = ^f[l + 3^ + 0((^f) 2 )]. (8) 
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Therefore the corresponding contribution to the beta function is 

Pi = -eg? + H9i) 2 H ■ (9) 

For the four-point function of type 3 , namely the one which involves a product 
of two Kronecker deltas of the form [r^]^ (5 a( j5f, c + $achd + dabbed), the situation is 
different since we now meet two kinds of diagrams, those with one A perloop and 
those in which the number of A's is equal to the number of loops plus one. On purely 
dimensional grounds we may decompose AT^ as 

Aif> (p) = 73 (p) + Hp) (io) 

with 

73 (p) = ^3-2 - + ••• 

e p e 

6 s (p) = + - (11) 

(i) Consider first the most infra-red singular contribution 5 3 {p). It comes from graphs 
with one A per loop,plus one. And its vertices can only be giS, since the oc- 
curence of one <?3 or more ( like the insertion of more A's) would entail some free 
replica summations vanishing with n. Furthermore the structure of those diagrams 
in [h{PiiP2,Pz,Pi)}ab is of the form 

[if } Wpi, • • • ,p 4 ) = (1 " <U)A 2 J d d q rS 4) ( Pl ,p 2 , q, Pl +p 2 - q) 
xG ab (q)G ab (p 1 +p 2 - q)T < i ) (q,p 1 +p2 - q,Ps,p 4 ) + permutations, (12) 

in which G a b(q) is the renormalized propagator. The integral over q is convergent, 
since (1 — 5 ab ) selects the part fo the propagator which falls off as 1/g 4 (up to loga- 
rithms). Therefore the renormalization of and of the propagators is sufficient to 
make the diagrams contributing to Ss{p) finite. No new renormalization is needed, 
and S 3 (pi,p2,P3,P4,) satisfies a Callan-Symanzik equation per se. 

(ii) 73 (p) comes from graphs with excatly one A per loop. They are all linear in g 3 
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since any higher power would again lead to free replica summations vanishing with 
n. Therefore the minimal renormalization of g%, defined as 

g 3 = ^g R [l + 2 9 -^ + 0((g R ) 2 )]. (13) 

is sufficient to make 73 (p) finite. The beta-function for g 3 follows : 

ft = -eg* + 2g R g R + 0(g R (g R f). (14) 

In general for the theory with the five coupling constants the same pattern governs 
the renormalization procedure. The functions and will both involve also 
terms with one A per loop, linear in g 2 and g%, which lead to a renormalization of 
g 2 and #3 ; they have both contributions with one more A which are made finite by 
the previous renormalization of . For one finds first terms linear in g 4 and 
quadratic in g 2 and g% which come from one A per loop ; they lead to a renormalization 
of g±. Then one finds terms proportional to A/p 2 , linear in g 2 and g$, and terms 
proportional to (A/p 2 ) 2 , which are made finite by the previous renormalization of 
Vi and of the propagator. For the situation is again similar except that there 
are now terms up to (A/p 2 ) 3 namely with k more A's than the number of loops, 
k = 0, • • • , 3. The k = terms lead to a renormalization of g§ ; the other ones 
are finite as a consequence of previous renormalizations. The five beta functions at 
one-loop order are then [|TJ : 

A = -m + 3gt 

A2 = -^92 + ^91(92 + 93) 

As = -eg 3 + ^91(92 + 9s) 

At = -e# 4 + 3<7i0 4 + 4(02 + 03) 2 

Ps = -^5 + ^94(92 + 93)- (15) 

The dimensional reduction fixed point, namely g\ = |e + 0(e 2 ) and g 2 = ■ ■ ■ = g$ = 
is unstable. (It is sufficient to notice that at this fixed point the matrix of derivatives 
1^- has an eigenvalue equal to |^ = — |e). 



4 Field renormalization 



At two-loop order a wave function renormalization appears. However a priori the 
random field introduces a privileged direction in the internal space along the unit 
vector v = ^(1, •••,!)■ We have thus to introduce the longitudinal and tranverse 
components of the fields 

-> 1 n 

(f) L = 0.v = J2<f>a (16) 

and 

— # — # 

0T = - 0L#- (17) 

A priori they are both renormalized by different factors and we define the renormal- 
ized fields through 

(f) L = \[Z~L$L, 4>T = \[Zt^T, (18) 

in which the n fields $ have finite correlation functions when e = (6 — d) goes to zero. 
In terms of these fields the quadratic terms in the action have the form 

PH = J ^{i[(V0 L ) 2 + (V0 T ) 2 ] - n^l}, (19) 

i.e. in terms of the renormalized fields 

(3H = J d d x{^\Z L {V^ L ) 2 + Z T {V^ T f}-nZ L ^ 2 L }. (20) 

This gives for the renormalized two-point function 

r$(p) = Z L (p 2 -nA)^ + Z T p 2 (5 ab - i) - S afe , (21) 

in which we have used the projectors along the vector v, the constant matrix whose 
elements are 1/n, and the projector transverse to v, the matrix 5 ab — 1/n ; E afe contains 
all the self-energy diagrams. 

From this expression one sees that it is expected that the wave function renor- 
malization of the longitudinal part of the propagator Z L renormalizes as well A, the 
variance of the random field. This is not obvious a priori; in the appendix an explicit 
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calculation to two-loop order is given, which shows that to this order Zt = Zl and 
that A does not acquire an independent renormalization. In a Langevin dynamics of 



non-renormalization of the temperature. 

5 Instability 

A derivation of the Landau theory with the five 4 coupling constants, was given 



distinct replicas, was attractive. We note that under renormalization g 3 is pushed 
further away from the origin. It is thus natural to examine whether this attractive 
interaction may lead to an instability. This question has been considered in some 
details by several authors |], |7|, §, [| . Here we wish to keep to very simple arguments. 

Therefore we now consider whether the Bethe-Salpeter kernel, for a pair of replicas 
of different indices, might develop a vanishing eigenvalue, thereby signaling bound 
state formation. Note that couplings other than g 3 vanish with n under iteration in 
the channel of two (distinct) replicas. The simplest iterative kernel (besides g 3 ) is a 
"bubble", i.e. a <5 3 -like contribution followed by two propagators, i.e. 



the same problem the non-renormalization of A is a natural consequence of the 



in [|TJ. It was shown there that the coupling constant g% = u^A, which couples two 




(22) 



An instability takes place whenever the spectrum of the operator 




(23) 



has a vanishing eigenvalue, in which 




(24) 



In position space this operator is 



h = [(V) 2 ] 2 - C~gl 



(25) 
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in which C is a positive constant. The four derivatives may balance the singularity 
at the origin of the attractive potential whenever 2d — 8 < 4. In that range the 
spectrum consists of bound states and positive energy scattering states. A simple 
scaling shows that the binding energies, in the domain d < 6 are proportional to 

4 _ 4 

| g 3 1 6 ~ d . Therefore whenever | g 3 \ 6 ~ d is larger than some critical value a bound 
state at zero energy appears, signaling an instability. Since gs is proportional to A, it 
means that, for A small, there is a small domain of size proportional to A 4// ( 6_d ) above 
the Curie line, in which the attraction between pairs of distinct replicas generate an 
instability. 

Quite generally it can be shown that the h operator of is, beyond the one- 
loop approximation, the Jacobian matrix of the Legendre transform || with respect 
to a source A ab (p) 

W(A ab ) + T(G ab ) = Wj2 A ab(p)G ab (p) (26) 

Z P ab 



i.e. 



h a b(p; p) 



d 2 Y 



dG ab (p)dG ab (p>) 

[Gaa{p)]- l [G bb {p)]- l 5 p+p , fl - rj(p, -p-p\ -p'). (27) 



On the other hand that double derivative has been identified with the replicon com- 
ponent of the Hessian matrix M, around the replica symmetric solution for the 
connected (G aa = G) and disconnected ( or connected through the random averaging, 
G a b — G) components. Thus the instability, i.e. negative eigenvalues of h, will force 
the emergence of a replica symmetry broken solution || . 

To conclude in one sentence, it may be said that the breakdown of dimensional 
reduction is not simply that the e-expansion needs to be corrected by non-analytic 
terms, but follows from a change in the phase diagram itself with the emeregence of 
a new, glassy, phase. 



Appendix : The propagator at two-loop order 
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The contributions to the self-energy due to the coupling constant g\ vanish at 
one-loop in the minimal subtraction scheme. At two-loop the singular terms of the 
self-energy are given by 

E ab = fi[A(p)6 ab + AB(p)] (A.l) 

with 

MP) = \l d d qi d d q 2 \- 1 (A.2) 
2 J q\q\ (p + qx + q 2 ) 2 

and 

B(p) = \l d d gi d d q 2 \- 1 (A.3) 
6 J qt <H VP + 5i + Q2Y 

9l = 9l (W/2)' (A ' 4) 
A lengthy, but standard, dimensional calculation in which d — 6 — e, leads to 

p 2+ 2e (1 _ e/2)2(1 _ e/4)2 r5(1 _ e /2)r(l + e) 



A = 



and 



£ = 



12e (l-3e/2)(l-3e/4) T(l - 3e/2) 
2 

~+ finite (A.5) 



P~ 2e (1 - 6/2) 2 (l - 6/4) 2 (l + 6) T 5 (l - £/2)r(l + 6) 

12e (l + e/2)(l + e/4)(l-e) T(l + e/2)r(l - e) 

= — + finite. (A.6) 

At this order the one-particle irreducible two-point function is given by 

r<J>(p) = z T (5 ab - -)p 2 + Z L -(p 2 - nA) - S a6 (A.7) 

and thus in the minimal subtraction scheme, we find that, since A and B diverge 
with the same ± ^ 

Z T = Z L = l + ^- + 0(gl), (A.8) 
and that indeed A does not require any additional renormalization. 
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